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SUMMARY 


A two-degree-of -freedom system, in which rolling oscillations are 
forced, is described. Details of the system, the theory, and the method 
of operation are discussed. It is shown that, although the system is 
characterized by nonlinear equations of motion, linearization of the 
equations by assuming small perturbations and constant coefficients yields 
sufficiently accurate results. 

The accuracy of the system was first investigated in tests under 
idealized conditions in which the aerodynamic derivatives were simulated 
by the action of gyroscopes and magnetic dampers. Later, the system was 
investigated in a high-speed wind tunnel using a simple model for which the 
aerodynamic derivatives could be estimated. 

The results of the tests showed that the quantities relating to the 
primary mode of operation of the system, that is, the rolling velocity 
derivatives, could be obtained satisfactorily. The same is true of the 
directional stability and the damping-in-yaw derivatives. Results obtained 
from the data-reduction equations for the rolling moment due to yawing 
velocity and due to sideslip angle were unreliable, however, and prevented 
the evaluation of these derivatives. 


INTRODUCTION 


The increasing emphasis placed on studies of the dynamic stability 
of aircraft has intensified the need for satisfactory methods for predict 
ing the dynamic stability derivatives. Most theoretical methods evolved 
to date are subject to limitations imposed to reduce the mathematical 
task to reasonable proportions. In order to obtain experimental checks 
of these methods and to provide sources of information in those areas in 
which theoretical methods are not yet available, there exists a need for 
experimental research on the stability derivatives. 

The experimental study of dynamic stability derivatives in wind tun- 
nels is, of course, not new. Many systems have been evolved for this 
purpose. References 1, 2 , and 3, for instance, discuss early 3ritish 
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systems in which the relative motion of the model and the air stream is 
generated hy model motion. A dynamic system recently developed at Ames 
Aeronautical Laboratory is described in reference 4; again the relative 
motion of the model and the air stream is generated by model motion. 

The rolling- and curved-flow wind tunnels now in use at the Langley Aero- 
nautical Laboratory of the WACA, for instance (refs. 5 and 6), are examples 
of other systems in which the model is held rigidly and the air flow gen- 
erates the relative steady yawing and rolling motion.. 

In the present case, a system in which the model executes forced 
oscillations with two degrees of freedom was developed. The details of 
its design and development are published herewith since they may prove 
useful to others who desire to develop dynamic systems. 


SYMBOLS AMD COEFFICIENTS 


The aerodynamic coefficients defined herein are referred to the 
system of body axes originating at the model center of gravity. The 
symbols and coefficients are defined as follows: 

b body length, ft 



f frequency, cps 
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Hx,Hy,H z 

moment of momentum about space-fixed axes 

hx > hy i hz 

moment of momentum about body axes 

lx 

moment of inertia about x axis, ft-lb sec 2 

I Z 

moment of inertia about z axis, ft-lb sec 2 

^xz 

product of inertia with respect to x and z axes, ft-lb sec 2 

I yz 

product of inertia with respect to y and z axes, ft-lb sec 2 

L 

rolling moment, ft-lb 

Lcp 

absolute value of spring constant of roll flexure pivots, 
ft-lb/radian 

Lq> 

rate of change of aerodynamic rolling moment with rolling 
angular velocity, ft-lb/radian/sec 

La 

9s 

equivalent viscous damping factor due to internal friction in 
roll flexure pivots, ft-lb/radian/sec 

LtJt 

rate of change of aerodynamic rolling moment with yawing 
angular velocity, ft-lb/radian/sec 

Lp 

rate of change of rolling moment with sideslip angle, 
ft-lb/radian 

L|3 

rate of change of rolling moment per unit rate of change of 
sideslip angle, ft-lb/radian/sec 

M 

pitching moment, ft-lb 


applied moments about space-fixed axes, ft-lb 

N 

yawing moment, ft-lb 

N$ 

rate of change of aerodynamic yawing moment with rolling 
angular velocity, ft-lb/radian/sec 

N\(r 

absolute value of spring constant of yaw flexure pivots, 
ft-lb/radian 


rate of change of aerodynamic yawing moment with yawing 
angular velocity, ft-lb/radian/sec 

N is 

equivalent viscous damping factor due to internal friction in 
yaw flexure pivots, ft-lb/radian/sec 

N(3 

rate of change of yawing moment with sideslip angle, 
ft-lb/radian 
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rate of change of yawing moment per unit rate of change of 
sideslip angle, ft-lb/radian/sec 

rolling angular velocity, radians/sec 

dynamic pressure, lb/sq ft, ^-PV 2 

maximum body radius, ft 

yawing angular velocity, radians/sec 

maximum cross-sectional area of body, sq ft 

roll input torque, ft -lb 

time, sec 

free-stream velocity, ft/sec 
space-fixed axes 

axes fixed with respect to the model (body axes) 
angle of attack, radians 

angle between free-stream velocity and longitudinal axis of 
support housing 

angle of sideslip, radians 

angle of roll, radians 

angle of yaw, radians 

phase angle between roll and yaw angles, positive for yaw 
leading roll, deg 

phase angle between roll angle and roll-input torque, positive 
for input ‘torque leading roll, deg 

angular frequency, 2itf, radians/sec 

mass density of air, slugs/cu ft 

amplitude of oscillation or out-of -phase component 
initial condition or in-phase component 

a( ) 

dt 
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ANALYSIS AND DESCRIPTION 
Coordinate System 


The coordinate system for which the equations of motion will be 
written is shown in figure 1. The x,y,z axes are fixed with respect 
to the model and are therefore body axes. The X,Y,Z axes are alined 
with the equilibrium position of the model and remain fixed in space. 

In order to illustrate the possible modes of motion of the model, 
a schematic drawing of the mechanical apparatus is shown in figure 2(a). 
As shown in this diagram, the main shaft (torque tube) of the apparatus 
is driven in a constant-amplitude oscillatory rotation about the X axis 
which is transmitted directly to the model. The axis of rotation of the 
torque tube is defined by the roll flexure pivots and ball bearing which 
support the tube within the support housing. Attachment of the model to 
the torque tube is through the yaw flexure pivots which allow the model 
a second mode of motion, that of angular rotation about the z axis. 


Equations of Motion 


In accordance with the Newtonian laws of motion, the time rate of 
change of the moment of momentum about axes fixed in space is equal to 
the summation of the externally applied moments about those axes: 


ZM X 


dHx 

dt 


dHy 

ZMv = — - 
1 dt 


ZM Z 


dH z 

dt 


/ 


( 1 ) 


Since the model is rigidly constrained in pitch, the second of equa- 
tions (l) need not be considered. 

The time rate of change of the moments of momentum about fixed- 
space axes in terms of the rate of change about the body axis system may 
be written (see, i.e., ref. 7)- 
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where 


3H X dh x 
dt - dt 


hy t 


d% 

dt 


dh 2 

dt 


hy q> 




( 2 ) 


h x - Cpl x + tl xz 

hy = -9Ixy “ Wyz ^ 

h z = cp^xz + 'flz 


(3) 


The applied moments and inertial reactions can then be expressed in terms 
of the body axis variables 


ZL - I x 9 + d xz ^ + I X y'l<'cp + Iy Z V 

ZN = Iz^ + Ixz9 ~ Ixy9 " dyz’P'l' 






Substituting the externally applied moments for the left-hand side of 
equations (4) gives 


T cos + 1 - 90,9 cos ^ + Lcp9 cos t 


- LtpCp + L4r v 4 + Lp [3 + Lp|3 
= I x 9 + Ixz^ + (Ixy9^ + Iyz’i' 2 ) 


(5) 




Ncp9 - + N3P + N p p = I z ^ + Ixz^ + ( " I xy 9 2 “ I yz c P'f ) 


A geometrical relationship exists between the angle of sideslip |3 
and the roll and yaw angles cp and t as follows: Let the angle between 

the free-stream velocity and the X axis be cti . 
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From sketch (a), 

w = V sin ^ Vi = V cos a ± 


and from sketch (b), 

v = w sin cp = V sin a x sin cp 



z 


From sketch (c). 


vi = v cos \|r = V sin ai sin cp cos \|r 


Sketch (a) 


y 2 es Vi sin \| r = V cos ai sin 
The angle of sideslip is then given by 





P = sin" 1 (sin a x sin cp cos \|r - cos ai sin \|r) 


\ 

\ 



The true angle of attack is also a 
function of the roll angle: 


a = 


sin 


-1 w 



= sin -1 (sin ai cos (p) 


( 7 ) 


Substitution of equation (6) and its 
time rate of change in equations (5) 
yields the equations of motion, and 
it can be seen that nonlinear terms 
will appear both in the aerodynamic 
moments and the inertial reactions. 



Sketch (c) 


In order that they be of practical value for use with experimental 
wind-tunnel data, linearization of the equations is necessary. If the 
approximation is made that 
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sin cp = cp sin ^ ^ 

cos Cp = cos i|/ = 1 

then the angle of attack becomes a = ctj., and the sideslip angle is given 

ty 

3 = cp sin a - \|r cos a ( 8 ) 

In general, the angle of attack is sufficiently small so that equation ( 8 ) 
may be written as 

3 = acp - t (9) 

and the time rate of change of sideslip angle is given by 

3 = acp - i (10) 


For large angles of attack for which these approximations are not valid, 
equation ( 8 ) and its time derivative may be used. 

If, in addition to the approximations made in equations (9) and (10), 
the second-order inertial terms in equations ( 5 ) are neglected, the equa- 
tions of motion may be written 


T + (L<p s + Lqj + 0 X 3)9 - (Lcp - 013)9 + (L\jr - LpH " Lp'k = *x9 + I xzt 
(Nq> + ccNj3 ) 4> + 0-1139 + (N^ s + Ity - Ng)| - (Ity + N^)\|r = I z f + Ixz? 1 


>( 11 ) 


The question immediately arises as to the accuracy with which equa- 
tions (11) represent the physical system. In order to answer this question, 
transient solutions of both the nonlinear and linearized equations of 
motion for a pulse -type roll disturbance T were obtained on the Reeves 
Electronic Analog Computer for conditions expected to be encountered in 
wind-tunnel operation. Comparison of the solutions showed that the only 
discernible effect of linearization was a negligible difference in the 
roll and yaw accelerations. From this result it was concluded that the 
use of equations (ll) would be a justifiable approximation, and that their 
use would yield data of sufficient accuracy. Further experimental veri- 
fication of the validity of these equations is given in a later section. 
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Solution of Equations of Motion 


In the following development the solution of the equations of motion 
will be given and arranged in a manner which will permit the determination 
of the desired aerodynamic parameters in terms of measured quantities. 

Yaw equation .- The second of equations (ll) may be rearranged and 
written 


I z * - ( w ts + ■ w |3)^ + (% + NpH = - x xzV + (Nq> + aWp )9 + ctNpCp (12) 


Equation ( 12 ) is a linear second-order differential equation, the solution 
for which is well known. The roll angle cp is forced to vary harmonically 
with time; 9= 9o sin wt. Then the solution, in general, consists of two 
parts: first, the complementary function, uniquely determined by the 

initial yaw angle and the initial yawing velocity \jr; and second, the 
particular integral, sometimes called the steady-state or frequency 
response, which is independent of the initial conditions, and is deter- 
mined only by the applied disturbance, the right-hand side of equation ( 12 ). 

The amplitude ratio and phase angle given by the particular integral 

are. 


*o 

% 


(a.Np + + (Ncp 


(% + N p - I z w 2 ) 2 + (% £ 


+ aWp) 2 w 2 
+ % - N p) 2 w 2 


(13) 


e 


tan -1 


(NA, + ctNo )(u 

- — + tan 

aNp + I xz io 2 


(N\jr s + Nf - Np)(o 

(N* + Np) - I z u 2 


Equation (l4) may be manipulated to give 


(14) 


aN 0 + Ixz 


0 )"- = 


w(Wq) + aNp) jjlty + Np - I z w 2 )cos e + (lity s + - Np)u sin ej 

(N\|f + Np - ) sin € - (N\|/ s + N\|/- - Ng)co cos € 


(Y9 i- al T ‘) - + Ixzt ° 2 ^0% + " IzW 2 )sin e - - Wp)w cos ej 

(% + Np - I z w 2 )cos e + (H^ s + fly - Np)u sin e 
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Substitution of these equations into equation (13) leads to two funda- 
mental relationships: 


\|r 0 ( N q) + 

*Po (K^ + Np - I z u 2 ) sin e - (N^ g + N^ - Np)w cos e 


(15) 


to aN P + I xz fa)2 

(Ny + Np - I z w 2 )cos e + (N^ s + N^ - Np)w sin e 


(16) 


The geometrical meaning of equations (15) and (l6) is apparent from the 
concept of the rotating vector diagram (see, e.g., ref. 8). Since 9 
is a harmonically varying function 9 = % sin <*>t, the right-hand side 
of equation (12) represents a harmonically varying yawing moment N 0 
with components proportional to the roll position, velocity, and accelera- 
tion. The vector diagram is shown in sketch (d) . The magnitude of the 

resulting yawing moment is 

1 *• No = % /(aNp + IxzW 2 ) 2 + (N<p + aNp) 2 w 2 

V », a/Vjyji*. 7 o'' 

Sketch (d) and angle £ is given by 


4 


(Nq, + aNn )w 

tan -1 1 

aWp + Ixz^ 2 


The particular integral shows that \|r is also a harmonic function, 
and equation (12) states that the sum of the moments R, resulting from 

the motion must be equal and 
opposite to N 0 * The vector diagram 
is shown in sketch (e). The magnitude 
of the resulting yawing moment N 0 
(which is equal to R) is 



N r 


= * 0 7(% + N p " I z a)2 ) 2 + ( N t s + Nyf - Np ) 2 w 2 


and the angle p is given by 


p = tan -1 


(N* s + N\|f - Np)to 
M + Np - I z w 2 
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The amplitude ratio and phase angle given by equations (13) and (ll) may 
then be written 


to = No/(Po 

^o No/to 


e = t + h 


Equation (12) may now be represented 
by combining the vector diagrams shown 
in sketches (d) and (e). It can be 
seen that equations (15) and ( 16 ) 
represent the summations of the out- 
of -phase and in-phase components of 
the vector diagram of sketch (f). 

The quantities in equations (15) and (l6) measured during operation 
of the dynamic apparatus are T 0 > to> u, and 6. At any discrete frequency 
these are insufficient data to solve the equations; however, use of the 
apparatus as a free -oscillation system provides additional data from the 
complementary solution of equation (12): 



^ 


\|r = S exp 


(% s + N* - Np) t 

. 2 T z 


♦i 


COS (x) a t - 


(Ky + N j - Np) 


2I z Wa 


sin (x) a t 


♦i 

77 — sin w Q t 


o>. 


“a' 


( 17 ) 


where the angular frequency of oscillation is given by 


w a = 


% + % / N t s + % - Np\ 

Iz V 2I Z ) 


( 18 ) 


For wind-off conditions, equations (17) and (l 8 ) may be written 


t = 



w n t 


2I z w n 


sin w n t 


+ — sin Un'l 
u n 


(19) 


where 



w n = 


( 20 ) 
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The contribution of the system damping to the oscillation frequency is 
negligible, and equation (20) may be written 


co 


n 



The moment of inertia may then be evaluated: 


Iz 


% 

4* 2 fn 2 


( 21 ) 


The system damping may be determined as 



(-) 

4'^t/ W i n d 0 ff 


(22) 


where t T and are ordinates of the envelope curve a time interval t 
apart; similarly, the aerodynamic damping becomes 


2I Z 

N* - Np = — 


2n 


V) 

^t/ w i n( i 


- 

on 


(23) 


and the directional stability parameter Np may be determined from 
equation (l8): 


Nr = 4* 2 I z f a 2 + 


(l^s + N(jf - Np)‘ 


4l, 


- N,, 


(24) 


Sufficient data are now available to solve equations (15) and (l6) for 
wind-on conditions : 


Ncp + ccNp 


+o 

9o 


(N^ + Np - I z co 2 )sin e - (N + N^r - Np)co cos e 


CO 


(25) 


Ixz 


^o (N\|f + Np - Izto 2 )cos g + (N\|f g + Nty - Np)w sin e 

cpo (J 2 


aNp 

C02 


( 26 ) 
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It should be noted that the product of inertia I xz may be obtained 

for wind-off conditions by elimination of the damping term in equa- 
tions (15) and (l6) to give 


+0 (% - I z u) 2 ) 
uftos € 


However, it was found that, for winged models at angle of attack, the 
slight elastic deformations of the model and support structure were suf- 
ficient to necessitate the use of equation (26) . 

Roll equation .- The first of equations (ll) may be rearranged and 
written as follows: 


1x9 “ (Lqfc + Lq, + aLg) 9 + (L(p - aL p )(p + I xz t - (I 4 - Lp)* + Lpt = T (27) 


Before proceeding with the total equation (27), we will consider the case 
of single -degree -of -freedom roll motion, defined by equation (27) when 
and its derivatives are identically zero. We proceed directly from the 
vector diagram given in sketch (g). 

Summation of the out-of -phase and in- 
phase vector components gives 

(Lq> s + hp + aLp)cp 0 w + T 0 sin A = 0 (28) 

Ix'fcw 2 - (L<p - aLp)9 0 + T q cos A = 0 (29) 

For wind-off conditions, equations (28) and (29) may be solved for 


0-4,+ 1-4+ <• 1-/3)+°“ 



T . 


— J 

'o Jjjf $0 Of* 


aLpjto I g 

Sketch (g) 


= 


T n sin A 


lx = 


cp 0 u 

T 0 cos A + Lcp9 0 

%w 2 


For wind-on conditions. 


Lq> + aLp = - 


T 0 sin A 


L P = a l L 9- 


cp 0 w 

T 0 cos A + I x 9oW' 

% 


9s 


(30) 

(31) 


(32) 

(33) 
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Equation (33) becomes indeterminate 



at zero angle of attack; however, 
another expression for Lp may he 
obtained from the more general equa- 
tion (27) which we now consider. 

The vector diagram representing 
equation (27) is shown in sketch (h). 
Summation of the out-of-phase and 
in-phase vector components gives 


T 0 sin A + I xz ^ 0 w a sin e + (Lcp g + L<p + aLp)cp 0 co + 


(L^ - L£)\|r 0 w cos e - Lp^ 0 sin e = 0 (34) 


T 0 cos A + (aLp - L<p)cp 0 + Ix^o^ 2 + I xz i|r 0 w2c0S e " 


(L^) - Lp ) \|r 0 co sin e - Lp^ 0 cos £ = 0 (35) 

Multiplying equation (34) by sin e, equation (35) by cos e, and adding 
and solving for Lp gives 


Lp = — 

acp Q cos g - to 

(Lcp s + Lcp + aLp)cpo w sin e 


( l^p — I x w ) cpo gos e — To cos(A — e) 


XZ T o' 


(36) 


where the product of inertia I xz is given by equation (26) . Multiplying 
equation (34) by cos €, equation (35) by sin e, and subtracting and 
solving for (L^ - Lp) gives 


H - Lp 



L<p + I x u 2 ) cp 0 sin e 


(L<p s + LA + aLp)cp 0 w cos e - 


T 0 sin(A - e) 


(37) 


Description of Apparatus 


A physical description of the drive apparatus and its various com- 
ponents is presented below. The next section deals with the functions of 
the instrumentation. 
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In discussing the drive apparatus, it is convenient to refer to the 
pictorial representation of the system presented in figure 2(h) . The 
model is allowed two angular degrees of freedom (rolling and yawing), is 
subjected to constant -amplitude forced oscillations in roll, and is free 
to respond in yaw. It is rigidly restrained in pitch (at a constant 
angle of attack from 0° to 15°) > and all externally applied pitching 
moments, lift, drag, and side forces are transmitted directly to the 
wind-tunnel structure. The apparatus is powered by a three-phase induc- 
tion motor mounted vertically over the support-strut housing at the rear 
of the apparatus. Rotary motion of the motor shaft is converted to an 
approximate harmonic motion of the main drive shaft (torque tube) by means 
of the eccentric block and drive flexure indicated in the figure. An 
enlarged drawing of these components is shown in sketch (i) to aid in 
visualizing the principle involved. 

Amplitude of the oscillatory rolling 
motion is determined by the amount of 
eccentricity in the eccentric block, and 
with the present system is 4°. The 
torque tube, which transmits the oscil- 
latory motion to the model, is supported 
within the stationary support sting by 
the roll flexure pivots at the forward 
end and by a ball bearing at the rear. 

Rigidly attached to the front end of the 
torque tube are the yaw flexure pivots 
on which the model is mounted. The only 
contact between the model and the appa- 
ratus is through the mounting at the for- 
ward end of the yaw flexures; sufficient 
clearance exists between model body and support sting to allow the model 
to yaw approximately 2° about the flexure axis. The maximum operating 
frequency of the system is approximately 20 cycles per second. 

To obtain accurate results with the system, it is necessary that 
rapid variations in oscillation frequency of the apparatus be eliminated. 
This purpose is served by the small flywheel shown mounted over the drive 
motor in figure 2(b). 

The use of flexure pivots in the apparatus resulted primarily from 
the need for avoiding the often undesirable damping characteristics of 
ball bearings. The flexures used were designed in accordance with infor- 
mation presented in reference 9> and were found to be entirely satisfactory 
from the standpoint of defining a rigid axis of rotation. 

Characteristics of flexure pivots .- The flexure pivots are designated 
as compensating type in reference 9 because, theoretically, the spring 
constant is essentially independent of load. The desirability of this 
feature in the present application is obvious, and the experimental veri- 
fication is shown in figure 3 for the case of the yaw flexure pivots sub- 
jected to lift loads of 0, 100, and 200 pounds. The invariance of the 
roll-flexure-pivot spring constant was equally good. 
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The internal damping in the yaw flexure pivots was obtained by means 
of the well-known oscillation-decay method. It was found that the rate 
of decay was not logarithmic but, instead, was a function of both ampli- 
tude and frequency. However, equation (22) was applied to the data in 
order to obtain an equivalent viscous damping factor N^ s . It was found 
that, as shown in figure k, the data could be represented by an equation 
of the form 


% 


s 


a + b^ 
f 


(38) 


Justification for this type of expression can be seen by converting to 
the form given in reference 10 where the internal damping is presented 
in terms of a "damping capacity" Aw, the damping work per cycle per unit 
volume. This quantity was found to be independent of frequency and to 
vary as the third power of the stress amplitude. Consider the case of 
constant-amplitude oscillation, t = sin cot. Then the damping capacity 
may be written (see, e.g., ref. r f) with the flexure-pivot cross-sectional 
area A and length h, 


Aw 


%s*o 2( ° 2 

Ah 



cos 2 o)t dt 


% s **o 2 a) 

Ah 


or, 


Aw = 


(a + bt 0 )2« 2 * 0 £ 
Ah 


(39) 


Although Aw was found to vary in reference 10 primarily as the third 
power of the stress amplitude (which is directly proportional to the 
oscillation amplitude) , examination of the results therein indicates that 
at the lower amplitudes the data may he approaching a second power vari- 
ation as indicated by equation (39)* Likewise, the damping capacity 
given by equation (39) is independent of frequency, and it is felt, there- 
fore, that the use of the equivalent viscous damping factor represented 
by equation ( 38 ) is justified. The data presented in figure 4 also show 

internal damping for the case of the yaw 
flexure pivots subjected to lift loads of 
0, 100, and 200 pounds. The application 
of an external load which produces no addi- 
tional damping is essentially impossible. 
However, it was found that the additional 
damping of a lift load simulated as shown 
in sketch (j) was negligible, and the data 
were obtained in this manner. Again it can 
be seen that there was essentially no effect 
of lift load. 



f 
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No successful method was devised for the application of a lift load 
which provided no additional damping about the roll axis. Consequently, 
the internal damping of the roll flexure pivots was assumed to be invariant 
with lift load, an assumption justified by the yaw-flexure results, and 
the fact that any small change that might occur would be a negligible 
percent of the aerodynamic roll damping. The equivalent viscous damping 
factor Lcp s for the roll flexure pivots was obtained by operating the 
apparatus under wind-off, tunnel-evacuated conditions, and it was found 
that these results could also be represented by an equation of the form 
of equation (38). 

The values of I<p s and N^ s were found to be approximately 20 percent 
and 12 percent* respectively, of the aerodynamic values obtained in the 
wind-tunnel investigation to be described later in the report. It should 
be noted, however, that the aerodynamic roll damping measured for the 
body-fin model used in the wind-tunnel test is much smaller than would be 
the case for a complete airplane model. 


Instrumentation 


Three quantities are measured by the system: angular roll and yaw 

positions of the model, and the roll -input torque applied to the model 
to maintain the forced oscillations in roll. These quantities are measured 
by means of resistance-type strain gages mounted on suitable beams, the 
locations of which are indicated in figure 2(b) . The amplitudes and 
phase angles of the oscillating strain-gage signals are obtained by pass- 
ing each of the signals through a sine-cosine resolver located at the 
drive motor and gear driven by the motor shaft (see fig. 2(b)) • The two 
output signals of each resolver are then directly proportional to the 
in-phase and out-of -phase components of the respective strain-gage signals 
referenced to the resolver rotor position. By the vectorial summations 
of these quantities the required amplitudes of oscillation are obtained, 
as well as the phase relationships necessary to establish the roll-input- 
torque and yaw-position phase angles with respect to the model roll posi- 
tion. Oscillation frequency is measured by means of an electronic counter 
driven by a portion of the oscillating roll-position strain-gage signal. 

Circuits . - It was noted in the preceding section that it is necessary 
to know the amplitudes of the roll angle, roll-input torque, and yaw angle, 
and the phase angles of the roll-input torque and yaw angle with respect 
to the roll angle. The desired quantities were measured by means of three 
identical electronic circuits, the essential features of which are shown 
in sketch (k) . 
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Sketch (k) 


The principle of operation of the circuits is as follows: A direct 

voltage is impressed on the strain-gage bridge which is to provide a 
measurement of the desired quantity. The resulting output is an alternat- 
ing voltage whose frequency is that at which the system is operating and 
whose magnitude is directly proportional to the amplitude of the desired 
quantity. This output is used to modulate a 400-cycle wave; the result 
is amplified and then passed through the resolver which reduces it to 
components in phase and out of phase with the position of the resolver 
rotor. These components are then demodulated and used to drive the gal- 
vanometers as shown in the sketch. 


Resolvers . - The resolvers utilized in the apparatus were of the 
standard induction type, commercially available. Two coils are accurately 
positioned at 90 ° with respect to each other on the stator and, with a 
1-to-l transformer ratio, produce two harmonic signals also phased 90 ° 
when the rotor is turned at a constant angular speed. 



The action of the resolver 
may be explained as follows: 

Assume the quantity to be measured, 
as represented by the oscillating 
strain-gage signal (hereafter 
referred to as the basic signal) , 
has an amplitude of oscillation, A, 
and phase angle, 5, with respect 
to the resolver rotor. This sig- 
nal is then defined as A sin(tot+5) . 
The resolver may be considered to 
provide two reference signals of 
unit amplitude, sin art, cos (*)t as 
shown in sketch (Z). Transmission 
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of the basic signal through the resolver effectively results in multipli- 
cation of the basic signal with the two reference signals: 


Thus the amplitudes of the resolver output signals are one half that 
of the basic signal, the frequencies are twice the frequency of operation 
of the system, and each of the resultant signals has been displaced, the 
first by an amount equal to one half the product of the basic signal and 
the cosine of the phase angle, the second by one half the product of the 
basic signal and the sine of the phase angle. 

Equations (4o) and (4l), which represent the components of the basic 
strain-gage signal in phase and out of phase with the resolver rotor 
position, also represent the time-varying voltages used to drive the 
galvanometers as shown in sketch (k) . 

Galvanometers . - A galvanometer element is essentially a single- 
degree-of-freedom system, and its amplitude and phase-angle responses to 
a harmonically varying forcing function are well known. If the frequency 
given in equations (40) and (4l) is much greater than the undamped natural 
frequency of the galvanometer element, the response of the element to that 
portion of the signal will be essentially zero. The first term on the 
right-hand side of each of equations (4o) and (4l) represents a voltage 
which is independent of time, however, and the galvanometer is capable 
of responding with a deflection G which is directly proportional to 
the static signal. The results, as shown in sketch (k), are static indi- 
cations of the in-phase and out-of -phase components of the harmonically 
varying quantity. 

Application of the proper calibration constants Co and Ci now per- 
mits the computation of the amplitude and phase angle of the desired 
quantity, the amplitude being equal to the square root of the sum of the 
squares of the in-phase and out-of -phase components, and the tangent of 
the phase angle being given by the ratio of the out-of -phase to in-phase 
components : 


Vi = A sin(cot + 5) sin cot 


= ^ cos 6 - — cos (2o)t + 5) 


(4o) 


Vq = A sin(o)t + 5) cos cot 


^ sin 5 + ~ sin(2cot -t 5) 


(41) 




( 42 ) 
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The phase angle thus obtained are those between the resolver rotors 
and each of the measured quantities; roll position; yaw position, and 
roll-input torque. It will be recalled, however, that the desired phase 
angles are those between the yaw position and roll positions, and between 
the roll-input torque and roll position. It is clear that these quantities 
are immediately obtainable from the above information if the relative phase 
angles of the three resolvers with respect to each other are known. Since 
it is convenient to the operation of the system to maintain these relative 
phase differences at small values, the resolvers are first alined with each 
other as accurately as possible. A calibration is then made from which the 
remaining phase angles are obtained. 


Effects of harmonics .- When the principle of operation of the resolver 
circuits was illustrated, it was assumed that the basic signal was a pure 
harmonic function. In general, the quantities to be measured are not pure 
harmonic functions, but may be distorted due to extraneous disturbances 
from support vibration and other sources, one of which was the slightly 
impure harmonic motion of the model produced by the eccentric block and 
drive flexure in converting the rotary motion of the motor into an oscil- 
latory motion of the torque, tube. In addition, the voltage output of the 
strain-gage bridge may have a mean value other than zero due to initial 
unbalance of the bridge. The voltage output may then be represented by 
the Fourier series 


00 

I (a * 


cos ncjt + b n sin ncjt) 


n=i 


or, in a more convenient form, 


where 


00 

Aq + ^ An sin(nojt + 6 n ) 

n=i 

A n = J a n 2 + b n 2 


c , -i b n 

6 n = tan — 


(43) 


In a similar manner, the resolver reference signals may not be pure har- 
monic functions; if, in addition, the two reference signals are not phased 
exactly 90 ° with respect to each other, but have some other phase angle 7 , 
the Fourier series representing the two signals may be written 

00 

^ B m sin mut (44) 

m=i 
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00 

^ B m sin (mart + 7) (45) 

m=i 


Then the in-phase and out-of -phase voltages used to drive the galvanometers 
will he given hy the products of equations ( 43) and ( 44) , and ( 43) and 
(45), respectively: 

00 

—T °° _«L 

Vi = A q ^ Bj^ sin mut + ^ A^B^ sin(rui)'b + ^n) sin mwfc 


m=i 


n=i m=i 


00 

= A 0 ^ B m sin mcjt + 

m=i n=i m=i 


1 1 ^f 5 {«» 


(n - m)wt + 6 n 


cos 


(n + m)o)t + 6 n 


(46) 


00 00 


V Q = A 0 ) B m sin(mcjt + 7) + ) ) A n B m sin(nwt + 6 n )sin(mcjt + 7) 


m=i 


n=i m=i 


Z oo 00 . _ _ 

/ \ v v AnBm r 

Bqj sin(mwt + 7) + ) ) — — i cos 

m=i 


n=i m=i 


cos 


(n + m)iot + &n + 7 


} 


(n - m)iot + 6>n ~ 7 

( 47 ) 


Examination of equations (46) and (47) reveals that time -independent 
voltages occur only when n = m, that is, for products of equal harmonics. 
Then the effective voltages, those to which the galvanometers will respond, 
may be written 


Vi 


Vo 


00 

I 

C ° S 2 5n A„Bn 

n=i 


00 

v 

cos(6 n - 7) 

L 

n=i 

2 


AnBn 


(48) 

( 49 ) 
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These equations show that the in-phase and out-of -phase components of 
each measured quantity may be in error by an amount equal to the sum of 
the products of equal harmonics contained in the quantity to be measured 
and in the resolver signals. If either of the signals has no harmonic 
content, but consists only of the fundamental (n = l), no error is caused 
irrespective of the distortion present in the other. Under normal wind- 
tunnel operating conditions it would be expected that, due to previously 
mentioned disturbances, some distortion would be present in the basic 
signals, and it is necessary that the resolver signals be free of all 
comparable harmonics in order to eliminate this possible source of error. 
For this reason the resolvers were driven by 90-tooth gears. A remaining 
possible source of low-frequency harmonic content in the resolver signals 
was slip ring "noise"; a visual examination of the resolver reference 
signals by means of an oscilloscope revealed no distortion, however. 

One other possible error due to the resolvers is indicated in equa- 
tion (49): the phasing of the resolver reference signals with respect 
to each other by some angle 7 other than 90°. For the resolvers used 
in the present investigation the angle 7 was accurate to 90°, ±l/2°, 
and the error from this source was, therefore, negligible. Now, since the 
fundamental amplitude, Bi, of the resolver reference signals is unity, 
equations (48) and (49) simplify to 


V ± = cos 6 X 


V n = — sin 8! 

The galvanometer deflections Gj_ and G 0 are directly related to the in- 
phase and out-of-phase voltages by the appropriate calibration constants 
Ci and Co* The amplitude and phase angle are given by 


Ai = 2 /co^Go* + Ci^Gi 


„ 2 , p . 2 


&i = tan 


.! C 0 G 0 

cIgT 


(50) 


Equations (50) are identical with equations (42) which were developed 
for the case of no distortion in the quantities to be measured. It can be 
seen that the system, in effect, filters out the unwanted harmonics and 
provides a measure only of the desired response. 
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Frequency characteristics ,- Mention was made of the alinement and 
calibration necessary to determine the final phase relationships of the 
resolvers. The method of alining the resolvers is as follows: While 

the oscillation frequency is maintained at some constant value, each of 
the three circuits is, in turn, connected to the output of the roll- 
position gage. In each case the appropriate resolver housing is then 
rotated until the out-of -phase component indicated on the galvanometer 
is approximately zero and is locked in this position. In- conjunction with 
these operations a check was also obtained on the frequency characteristics 
of the measuring circuits. It should be recalled that the phase angles 
necessary for data reduction are measured with respect to the roll posi- 
tion. Then it is clear that the absolute value and frequency dependence 
of the phase angles indicated by the resolver circuits is immaterial if 
these characteristics are identical for each circuit. The procedure was 
as follows: After alining the resolvers, each of the circuits was, in 

turn, connected to the output of the roll-position gage, and the ampli- 
tude and phase angle were computed for discrete frequencies covering the 
available range of operation. It was found that the phase shift varied 
slightly for the three circuits, the phase angle differences varying in 
a manner linear with frequency. The phase shift difference was then used 
as a calibration for correcting the phase angles obtained during testing. 


Operation of the System 


The operational procedures followed in conducting wind-tunnel inves- 
tigations with the testing apparatus are described briefly below. It 
should be noted that the procedures are essentially identical when the 
apparatus is used for calibration tests such as are described in a later 
section. 

(1) With the model rigidly restrained in yaw, measurements of roll 
position cp, roll-input torque T, and the roll-input-torque phase angle 
A are obtained for single-degree-of -freedom roll oscillations, at dis- 
crete frequencies throughout the available range, by means of the resolver 
circuits previously described. Data thus obtained, for both wind-off 
tunnel-evacuated, and wind-on conditions, and the calibrated roll-flexure 
spring constant, L<p, provide Ly Q and I x from equations (30) and (3l)> 

L^p + aLp from equation (32), and, for angle of attack other than zero, 

Lp from equation (33)* 

(2) With the model free to yaw, the model is driven in roll at the 
yaw natural frequency. When the model has attained sufficient yawing 
amplitude, the drive motor providing the oscillatory rolling motion is 
stopped abruptly, resulting in a single-degree-of -freedom free-yaw oscil- 
lation. A time history of the decaying yawing motion is obtained by 
means of a recording oscillograph. Again these data are obtained for 
both wind-off tunnel-evacuated and wind-on conditions, and with them and 
the calibrated yaw-flexure spring constant, N^, equations (21) through 
(2k) yield I z , N^ g , - Np, and Np, respectively. 
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(3) With the model free to yaw, the apparatus is operated as a two- 
degree -of -freedom system at discrete frequencies throughout the operating 
range. At each frequency measurements of the roll position cp, yaw posi- 
tion \|r, roll-input torque T, and the phase angles e and A are obtained 
by means of the resolver circuits. Use of these data in conjunction with 
those obtained from (l) and (2) above provides + aNp and I xz from 
equations (25) and (26), Lp and - Lp from equations (36) and (37)* 


EVALUATION OF SYSTEM 


It should be noted that the aerodynamic moments produced by angular 
velocities are, in general, much smaller than the moments produced by 
angular displacements, and the disturbances inevitably present in any 
wind tunnel may be of the order of magnitude of the quantities to be 
measured. An analysis of the manner in which the effects of these dis- 
turbances are eliminated in the present system has been discussed in the 
section on instrumentation. The results of such analyses are not suf- 
ficient information to determine completely the over -all accuracy; however, 
a reliable indication of the accuracy can be obtained where it is possible 
to simulate the quantities to be measured with known inputs to the system. 
The results of tests of this type will be discussed. A further check of 
the system accuracy when operating under wind-tunnel conditions was obtained 
through the testing of a simple configuration for which, it was felt, cer- 
tain of the aerodynamic derivatives could be calculated fairly accurately. 
The results of these tests are discussed at the end of this section. 


Evaluation of the System With Simulated Aerodynamic Derivatives 


Cross derivatives N<p and L\j).- Tests with the simulated aerodynamic 
derivatives were made in order to assess the accuracy of the system and 
the validity of the assumptions made in the linearization of the equations 
of motion. The aerodynamic cross derivatives N(p and were simulated 
by means of a variable-speed gyroscopic motor as indicated in sketch (m) . 

The torque due to precession of a gyroscope is 
well known (see, e.g., ref. 8), and it is this 
property which was utilized in the present appli- 
cation. The gyroscope was mounted in place of 
the model, and, with the rotor oriented as shown, 
the simulated cross derivatives provided by the 
gyro are given by 


L$ = -Ncp = IQ 



Sketch (m) 
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where 

I rotor moment of inertia 

Q angular velocity of rotor 

The results of the gyroscope calibration are shown in figure 5 where 
the range of values given is expected to cover all possible wind-tunnel 
results. It can be seen from figure 5 that the accuracy with which Nq> 
is determined is entirely satisfactory. It was found that the values 
of L\Jr obtained from equation ( 37 ) were unreliable, however, and these 
results are not presented. This was also found to be true for the calcu- 
lation of Lg from equation ( 36 ) and a re-examination of the method of 
determining these derivatives will be necessary. 


Damping-in-roll derivative Lm.- The results of tests in which the 

was simulated are shown in figure 6 . This 
indicated in sketch (n). 


aerodynamic roll damping I<p 
derivative was simulated in the manner 
proportional to the angular velocity 
<p was produced about the x axis 
by means of a magnetic damper acting 
through a flexure-pivot universal 
joint on an arm connected rigidly to 
the roll torque tube. The damping 
moment was generated by operating 
the system with a direct-current 
field excitation and a short-circuited 
armature. If the current generated 
in the armature and the armature 
resistance are designated as I and R, 
respectively, the simulated aerodynamic 
derivative may be expressed, for small 
angular displacement (p, as 


A moment 


L<p = 


2jtcI 2 R 
^22 
9o W 



where c is the conversion from joules to inch-pounds. With the damper 
inoperative, an equivalent damping factor was obtained for the internal 
friction in the roll flexure pivots, universal joint, and flexible mount- 
ing of the damper armature. This quantity was subtracted from the damper - 
on data to give the results shown in figure 6 . These results show the 
accuracy with which the system measures the simulated aerodynamic deriv- 
ative L(p in the absence of tunnel disturbances. 


Damping in yaw and directional stability derivatives 


Cn r ~ ^nfi and 


C n p • - No satisfactory method was found for introducing simulated values 

of aerodynamic yaw damping and directional stability into the system 
during the previously described tests, and it is, therefore, necessary 
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to analyze the results of the wind-tunnel tests, described in the follow- 
ing sections, to determine the accuracy with which these derivatives were 
obtained. The procedures for obtaining and reducing free-oscillation data 
are well known, as are the factors affecting the accuracy of these data. 

The uncertainties in each of the measured quantities present in the equa- 
tion defining Cn r - Cn3 and Cnp have been calculated. The square root 
of the sum of the squares of these uncertainties indicated over-all 
uncertainties of the order of ±0.055 and ±0.057 for single values of 
Cn r - Cno and Cnp, respectively. Not accounted for in these values, 
however, are possible effects of random air-stream disturbances and support 
vibrations which would appear as scatter in the data. The standard devia- 
tion of a number of observations at a given Mach number was, therefore, 
calculated; this was found to be of the order of ±0.127 for Cn r - Cn^ • 
and ±0.028 for Cnp* 


Evaluation of the System in the Wind Tunnel 


The tests described previously, in which the system accuracy was 
assessed from the standpoint of ability to measure simulated aerodynamic 
moments, were necessarily conducted under somewhat idealized conditions. 

It is, of course, impossible to simulate all of the loading conditions 
and disturbances to which the model is subjected in the wind tunnel, and 
these tests provide a measure of the accuracy of the system when there 
are essentially no extraneous disturbances present. It has been shown 
in the analysis of the previous section that these disturbances, theoret- 
ically, cause no error; however, it is desirable to check the capabilities 
of the system operating in the wind-tunnel environment for which it was 
designed. 

An evaluation of the system in the wind tunnel necessarily depends 
on testing a model for which the aerodynamic derivatives may be calculated 
accurately. It is, of course, realized that an exact check on the system 
accuracy is not possible by this means since present theoretical knowledge 
is insufficient to permit exact calculations of the lateral dynamic sta- 
bility derivatives of fin-body combinations; a serious decrease of the 
system accuracy due to tunnel disturbances should, however, be evident. 

For the tests described below, a simple model was constructed for 
which, it was felt, the rolling-velocity derivatives Cj p and C np could 
be estimated with a fair degree of accuracy. The model, which consisted 
of a vertical triangular fin mounted on a cylindrical body with an ogival 
nose (see fig. l) , was tested with the dynamic system in the Ames 6- by 
6-foot supersonic wind tunnel over a range of subsonic and supersonic 
Mach numbers at zero angle of attack. 

Estimation of stability derivatives .- Estimation of the stability 
derivatives was based on existing theoretical methods adapted to the 
particular model used in the wind-tunnel tests. The following is a brief 
description of the methods by which the derivatives were estimated. 
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In determining the roll-velocity derivatives, and C n p, it was 

assumed that the fin acts as a portion of a rolling triangular wing as 
indicated in sketch (o). The solution for the resulting pressure dis- 
tribution is well known (see, e.g., ref. 11) and was used directly to 
calculate the damping-in-roll derivative and the contribution of 


the fin to the cross derivative C n p 
loading on the fin of course provide 
no rolling moment, but contribute 
directly to the yawing moment, and, 
therefore, must be considered in 
estimating C np tor the fin-body 
combination. Methods are available 
for computing the interference load 
distribution for fin-body combinations 
at supersonic speeds (refs. 12, 13 , 
and l4) but these are, in general, 
very laborious and for this reason 
the method described below was applied 
in the present case. 


Pressures induced on the body by the 



\ 


Sketch (o) 


It was assumed that the pressures induced on the body could be cal- 
culated as though the body were a flat plate at zero angle of attack 
adjacent to the fin. Reference 15 
gives the solution for a case similar 
to this, namely the pressure distribu- 
tion induced on a flat plate at zero 
angle of attack adjacent to a tri- 
angular sector at a constant angle 
of attack. The angle-of -attack dis- 
tribution on the fin may be repre- 
sented as the sum of a constant angle 
Pr N 

increases linearly in the spanwise 
P z ^ 

a = — j as shown m 


of attack a = — and one which 


direction 
sketch (pK 


The pressure induced on the 
body by the constant angle -of -attack 
distribution is given directly by the 
solution in reference 15* The induced 
pressure resulting from the effective 
twist due to rolling velocity may be 
found by integrating across the span 
of the wing. This integration cor- 
responds to the superposition of an ^ ^ 

infinite number of the lifting sectors Sketch (p) 

along the span, each sector having an infinitesimal angle of attack. The 
yawing moment due to the induced pressures on the body was added to that 
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due to the loading on the fin to obtain the total yawing moment due to 
roll velocity. Because of the approximations necessary in the above 
estimate, the yawing moment due to rolling velocity was also calculated 
using the slender-body theory of reference l6 in order to provide a 
check of the order of magnitude. 

In estimating the yawing derivatives it was assumed (a) that, in 
lieu of a more exact theory, the yawing moment of the body is that given 
by slender-body theory which, of course, is known to be incorrect at the 
higher supersonic Mach numbers and (b) that the yawing moment of the fin 
due to yawing velocity is that given by the moment, about the appropriate 
axis, of the loading on a portion of a pitching wing. Sketch (q) shows 
the portion of the wing referred to as a cross-hatched area and the 

pitching axis as the axis of yaw. The load- 
ing was calculated according to reference 17 
for the arrangement as shown to determine 
the yawing moment due to yawing velocity. A 
similar method was used for estimating the 
yawing moment due to sideslip angle. The 
slender-body-theory results were obtained 
from reference 18. The above methods, of 
course, applied only to the derivatives for 
supersonic Mach numbers. No methods were 
available of sufficient accuracy to permit 
the estimation of the derivatives at subsonic 
Mach numbers. 

Comparison of estimated stability derivatives with measured values . - 
During the wind-tunnel tests the rolling derivatives C^p and C n p were 

obtained over a range of frequencies from 7 to 20 cycles per second. No 
significant variation of these quantities with frequency was found, and 
the values which are presented in figures 8 and 9 are the arithmetic 
mean of the values obtained. The values of Cn r - Cnp and C n p are those 
which were obtained at a single frequency through free-oscillation tests. 

A comparison of the estimated and measured rolling derivatives is given 
in figures 8 and 9 tor the supersonic Mach numbers investigated. The 
measured values are also shown for the subsonic Mach numbers. The agree- 
ment between the experimental and estimated values indicates that there 
were no serious effects of wind-tunnel disturbances on the system accuracy. 
Figure 10 shows the measured and estimated damping in yaw. Here the 
agreement is poor, particularly at the highest Mach number. The same is 
true of the directional stability parameter Cnp.? shown in figure 11, 
where the agreement is satisfactory except at tne highest Mach number. 

The methods used in estimating the theoretical values of the derivatives 
may not be reliable and have not been adequately checked experimentally. 

A large part of the discrepancy between theory and experiment is believed 
to be due to deficiencies in the theory. 
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CONCLUSIONS 


A study of a system designed to measure the static and dynamic 
lateral stability derivatives has lead to the following conclusions: 

1. Tests in which aerodynamic derivatives were simulated have shown 
that the apparatus correctly measures known inputs of the order of magni- 
tude anticipated for airplane models* 

2. The investigation has shown that the aerodynamic derivatives 
associated with the primary mode of operation of the system, that is the 
rolling derivatives, and the directional-stability and damping -in -yaw 
parameters may be obtained with satisfactory accuracy. The data-reduction 
equations from which the rolling moment due to yawing velocity and due 

to sideslip angle were determined yielded results which were unreliable 
and prevented evaluation of - Ci^ and C^. 

3. Tests of the apparatus utilizing a simple model showed that the 
wind-tunnel environment apparently had no deleterious effects on the 
system accuracy. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 24, 1954- 
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Illustration of coordinate system. Displacements and moments are shown in the positive 
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(a) Schematic drawing* 

Figure 2.- Illustration of dynamic testing apparatus. 
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Figure 2.- Concluded. 
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Figure 3»- Effect of simulated lift load on yaw flexure pivot spring constant. 
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Figure 4.-* Effect of simulated lift load on the internal damping of the yaw flexure pivots. 
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Figure 5.- Comparison of gyroscope input and measured output of the simulated aerodynamic cross 

der i vat i ve l^p . 
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Figure 6.- Comparison of magnetic -damper input and measured output of the simulated aerodynamic 

damping-in-roll derivative L<p. 
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Figure J.- Model used for wind-tunnel tests. 
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Figure 8.- Variation with Mach number of experimental and estimated values of damping-in-roll 

coefficient, C i . 
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Figure 9»- Variation with Mach number of experimental and estimated values of cross -derivative 

coefficient/ C n p« 
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Figure 10.- Variation with Mach number of experimental and estimated values of damping -in -yaw 

coefficient, C nr - C n £. 
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Figure 11.- Variation with Mach number of experimental and estimated values of directional stability 

coefficient, C n p. 
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